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Abstract: 
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possible anyon-like excitations in the spectrum. It is also intimately connected to the breaking 
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1 Introduction 



Chiral symmetry and its realization at finite temperature have been matter of intense studies 
for a long time. One of the most extensively used laboratories for this purpose has been the 
Gross-Neveu model [1] both at zero and finite temperature [2] [3]. In the large- N limit the 2 + 1 
dimensional version of this theory is renormalizable and provides a convenient arena for the 
survey of chiral symmetry; for instance, its dynamical breakdown at zero temperature and its 
high temperature restoration. Recently, there has been also considerable interest in studying 
the properties of chiral symmetry at finite density, or equivalently in the presence of a chemical 
potential [4-6]. This is believed to be relevant to the understanding of the physics of hot and 
dense matter which is expected to be probed in the laboratory by the RHIC experiments [7]. 

In this work, we address the thermodynamics and chiral properties of the 2+1 dimensional 
Gross-Neveu model in the canonical formalism. The use of the canonical formalism follows 
naturally when one introduces a constraint on the fermion number. Our discussion, however, 
will be based on a generalized representation of the canonical partition function that encodes 
the presence of a real chemical potential as a special case. It is shown that the approach 
yields similar results to the grand canonical formalism outcomes, the latter being the natural 
framework for studying systems in the presence of a chemical potential [8,9]. This generalized 
partition function will also be studied for imaginary values of the chemical potential where it 
corresponds, rather surprisingly, to a real free-energy density. 5 

We argue that the presence of an imaginary chemical potential is intimately connected to 

the possibility of having anyon-like excitations in the spectrum of the theory. We also show that 

the imaginary chemical potential emerges naturally in the 2+1 dimensional Gross-Neveu U(l) 

gauged model at finite temperature. Such a theory has infinitely many Z-vacua [13] around 

which the gauge field may fluctuate. These fluctuations have been shown to be connected to 

peculiar excitations whose charge is not an integer multiple of the elementary charge of the 

theory. We interpret these excitations as being anyon-like [14] and argue that our generalized 

canonical partition function gives their free-energy density. As the gauge field fluctuations 

around the Z-vacua lead to the spontaneous breakdown of the discrete Z-symmetry, it appears 

that the latter is tied to an any on confinement/ deconfinement transition. Finally, we study the 

chiral properties of the Gross-Neveu model in the presence of an imaginary chemical potential. 

In this case the theory is chirally symmetric at high enough temperatures but it appears to be 

unstable at T = 0. For any non-zero value of the imaginary chemical potential, chiral symmetry 

is broken at a certain temperature as the system cools down. This corresponds to a second-order 

phase transition where the free-energy density of the system is given by a remarkably simple 

analytic expression. There also exists a particular value of the temperature at which the system 

5 The study of the thermodynamics of fermionic systems in the presence of an imaginary chemical potential 
has been considered in the past [10, 12, 11]. Such studies, however, have been primarily focused on the technical 
virtues of the imaginary chemical potential formalism in lattice simulations. 
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appears to be in a bosonic phase. 

The paper is organized as follows. In Sec. 2 we briefly review the 2+1 dimensional Gross- 
Neveu model at finite temperature. In Sec. 3 we set the stage for studying the model with a 
constraint on the fermion number. Using a generalization of the well-known formula which gives 
the canonical partition function [10, 12, 11] we derive analytic expressions for the free-energy 
density, the fermion number density and the chiral order parameter. Our formalism coincides 
with the standard grand-canonical approach when the chemical potential takes real values. In 
Sec. 4 we study the model in the presence of an imaginary chemical potential and show that it 
corresponds to a system possessing anyon-like excitations. We establish a connection with the 
breakdown of the Z-symmetry [13] of the U(l) gauged Gross-Neveu model at finite temperature. 
In this way we give a physical interpretation to the canonical partition function with imaginary 
chemical potential as representing the free-energy needed to immerse an excitation of imaginary 
charge in the spectrum. We argue that such excitations are anyon-like. Next, we present results 
related to the properties of chiral symmetry in the presence of an imaginary chemical potential 
at finite temperature. In Sec. 5 we summarize and discuss possible implications of our results. 



2 The 2+1 dimensional Gross-Neveu model at finite temperature 



The Euclidean Gross-Neveu model may be defined by the Lagrangian density 6 [2, 15] 

C = -$ i frl> l -±(r'c 1 )-. (1) 



where if) 1 , ip % , i = 1,2, .., N, are four-component Dirac fermions and g is the coupling. In d = 3 
the massless model above is [7(27V)-invariant and possesses a Z2 discrete "chiral" symmetry 



(2) 



This model has been extensively used as a testing ground for studying the mechanism of chiral 
symmetry breaking in QCD [6]. 

For large- N, the model is studied in a 1/N expansion where it is renormalizable [15]. One 
introduces an auxiliary scalar field \(x) and integrates the fermions out. The partition function 
(generating functional) reads 



Z = 



J (VX)e- NI eff^ , 



(3) 



/ \ 

a' 



with a', i = 



6 For the Euclidean gamma matrices we use the Hermitean representation ji — 

\° ° J 

1,2,3 the usual Pauli matrices and 70 = 73. In this reducible representation there exist two 75-like matrices 



1 



and 



/ \ 

1 



-1 
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I ef f(X,G) = ±Jd 3 x\ 2 (x)-Tr[H$ + \)}, 



(4) 



where the rescaled coupling G = gN is kept finite as N — ► oo. 

To leading order in l/N the renormalized theory manifests itself in two different phases 
distinguished by a zero or a non-zero expectation value of (A). This manifestation depends 
on the value of the renormalized coupling l/G r as compared to the critical coupling 1/G* = 
4(2-7r)~ 3 / d 3 p/p 2 . For 1/G r > 1/G*, (A) = and the theory is in a weakly-coupled phase 
where chiral symmetry is unbroken. For l/G r < 1/G*, (A) ^ and the theory is in a strongly- 
coupled phase where chiral symmetry is broken. In the latter case the fermions acquire a mass 
proportional to (A). Clearly, (A) plays the role of an order parameter for the chiral phase 
transition. 

From the Euclidean formulation above one straightforwardly switches over to thermodynam- 
ics at temperature T = 1/(3 by making xq finite with length L = 1/T 7 and imposing periodic 
(antiperiodic) boundary conditions over the interval [0, L] for bosonic (fermionic) variables. In 
this way the bulk theory in 2+1 dimensions corresponds to a two-dimensional quantum system 
at zero temperature. It then follows that the T = system can be "prepared" , by appropriately 
tuning the coupling constant, to be either in the chirally symmetric or in the chirally broken 
phase. Had it been "prepared" to be in the broken (ordered) phase, one would expect that 
there exists a high temperature phase transition to the symmetric (disordered) phase. Such a 
transition is allowed (not forbidden) by the Mermin- Wagner- Coleman theorem as the relevant 
symmetry is discrete. 8 

To study the chiral symmetry restoration one calculates the partition function (3) by the 
steepest descent method for large- N. This amounts to performing a l/N expansion for Z around 
its saddle-points, which correspond to uniform values of (A). The latter are obtained from the 
gap equation 



This expression is divergent but all UV divergences at finite temperature are the same as at zero 
temperature [15]. Therefore, one can renormalize (5) by substituting for 1/G its corresponding 
renormalized value at T = 0. For the system to be in the chirally broken phase at T = one 
sets 



7 In the following we shall invariably use T, (3 and L bearing in mind the relations (3 = L = 1/T. 

8 Compare this situation to the 2+1 dimensional O(N) vector model [16, 17]. In that case, the bulk theory 
again corresponds to a two-dimensional quantum system at T — which can be in either an 0(iV)-symmetric 
phase or an 0(iV)-broken phase. However, the T > theory can only be in the 0(iV)-symmetric phase as the 
relevant O(N) symmetry is continuous and cannot be broken in d = 2. 




(5) 
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where m > is the dynamically induced mass of the elementary fermions at T = 0. Then, from 
(5) one obtains 



<A) = I ln 



e Lm/2 _|_ ^/gLm _ 4 



(7) 



which gives the dependence of the order parameter (A) on the temperature. In particular, (A) 
vanishes at the second-order phase transition point 1/L C = T c = m/21n2 where chiral symmetry 
is restored. At this point and for higher temperatures, the free-energy density is given to leading- 
N by 9 

/ _ 303) 
N 2vrL 3 ' 

and coincides with the free-energy density of 2N massless four-component Dirac fermions [19]. 



(8) 



3 Thermodynamics of the 2+1 dimensional Gross-Neveu model in the canon- 
ical formalism 



3.1 General Setting 

The canonical formalism for the analysis of the thermodynamics of a system has been recently 
employed in studies of fermionic systems at finite baryon density [11, 12]. The reason is that 
it bypasses the usual sign problem of the Euclidean fermion determinant for real values of 
the chemical potential. Since the three-dimensional Gross-Neveu model does not undergo this 
sign problem, it has been extensively studied in the standard grand-canonical formalism both 
analytically and numerically [4, 5]. Nevertheless, it would still be interesting to perform a direct 
analysis in the canonical formalism in order to test results obtained previously. Furthermore, in 
doing so, we obtain some surprising new results which shed new light into the thermodynamics 
of the model. 

The canonical partition function can be obtained as the thermal average over eigenstates of 
the number operator N = J d 2 x t/^(t, x)?/>(t, x) with fixed eigenvalues B. Namely, [8] 

Z(T,V,B) =Tr(e- /3A 5{N-B)^ , (9) 

where H is the Hermitian Hamiltonian. As usual, N measures the excess of fermions over 
anti-fermions in the spectrum. 

If one anticipates that for certain physical conditions the spectrum contains free fermions or 
anti-fermions, then B must be an integer. This then leads to the following representation for 
the canonical partition function [10] 

Z(T, V, B) = — / dcf> e-^ B Z(T, V, i<f>) , (10) 

27T JO 

9 In our calculations we always normalize the free-energy density in such a way that it vanishes in the bulk, i.e. 
at T=0 [18]. 
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where, Z(T, V, icp) is the grand-canonical counterpart with imaginary chemical potential. If on 
the other hand one is interested in physical situations where B is the mean fermion (anti-fermion) 
density, then B is real, not necessarily integer, and (10) is no more a valid representation. 

In this work we propose a generalization of Eq.(10) which is suitable for both analytic 
and numerical work. Our representation includes both cases of integer and non-integer real 
values of B. Furthermore, it leads to some unexpected new results when B is imaginary. Our 
representation follows from Eq.(9) if we write the delta- function constraint with the help of an 
auxiliary Lagrange multiplier scalar field 9{t) as 



Z(T,V,B) 



(VX)(V6)e~ NX ^f (x ' 9 ' 9 ' B) 

■0F(T,V,B) 



(11) 



which enforces the averaged fermion number constraint (N) = B. Here, F(T, V, B) is the free- 
energy and X e f f is the following effective action 

T e// (A, g; 9, B) = i^ £ dr 9(r) + ^ £ dr J d 2 x A 2 (r, x) - Tr[ln(p + A(r, x) + i 7 o0(r))] L • 

(12) 

To evaluate the free-energy in (11) for large-iV, we expand I e ff around its stationary points 
assuming constant (translation invariant) (A) and (9) configurations. These satisfy the following 
set of saddle-point equations 10 

d 2 p 1 



1 

G 

B 



E 



(27r) 2 p2 + K + (0)) 2 + (A) 2 ' 



lim — 

p^o L 



d 2 p 



E 



p2 + K + <#»2 + (A>2' 



(13) 

(14) 



where oo n = (2n + l)vr/L, n = 0, ±1,±2, ... [20,19]. We can subtract the UV divergences in 
Eq.(13) [19] by adjusting the renormalized coupling as in (6). At the same time the system at 
T = can be arranged so as to break chiral symmetry. After some algebra we obtain 





B 

N 



<A>- 
2vrL 



1 

m + — 

1j 



ln( 



1 + e 



-L(\)-iL(0)\ 



+ 



ln( 



1 + e 



-L(X)+iL(e)\ 



In 



1 + e 



-L{\)+iL(8) \ 



1 + e -L(A>- 



TU8))- Ll2 (" 



Li 2 (- 



-L(\)-iL(ey 



(15) 
(16) 



where B = B/V is the fermion density. From (11) and (12) we can now calculate the renormal- 
ized free-energy density / = F/V to leading- N and obtain 

(A) 3 m(A) 2 



!/(/?, m; (A), B) = _ 



+ 2iB{0) + 



1 



2vr " 6vrL 3 
ln(e- L W) (Li 2 (-e- L W+ iL W) + Lh{-e' L ^- iL ^ 



Li 3 (_ e -£(A)+iiW) + Li 3 (- e - L < A >- iL W) 



(17) 



10 As discussed in [19, 20] the regulating parameter p in (14) is necessary in order to take care of the fact that 
the Tr[...] and ^[...] operations do not commute. 
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The functions Li n (z) are the standard polylogarithms [21]. 

In principle, the free-energy density (17) together with the gap equations (15) and (16) are 
sufficient for studying the thermodynamic properties of the Gross-Neveu model to leading- N. 
Before we proceed, however, we observe that requiring the free-energy density in (17) to be real 
we are naturally led to distinguish two cases. Namely, as the terms involving polylogarithms 
in (17) are real for both real and imaginary values of (9), we can have either 1) (9) = imaginary 
and f?=real, or 2) (0)=real and S=imaginary. Clearly, 1) would correspond to the usual case 
of the Gross-Neveu model with real chemical potential. Nevertheless, we will show that case 2) 
probes some interesting properties of the Gross-Neveu model too. 



3.2 Relation to the grand-canonical formalism 



For (#)=imaginary (17) is the free-energy density of the Gross-Neveu model with ordinary 
(real) chemical potential. Indeed, setting i(9) = /j, < m we find that (15) and (16) coincide 
with the corresponding expressions for the gap equation and the fermion density, the latter was 
obtained for the first time in [5]. Namely, 



(A) = 

B 

N 



2vrL 



2cosh(L/i) + J(e Lm - 2cosh(L / u)) 2 - 4 



In 



1 _|_ e -L(\)+Ln" 
I _|_ e -L(\)-Lfj. 



Li 2 (- 



+ 



Lh (- 



-L(\)-L^ 



(18) 

(19) 



The thermodynamic properties here are well-known. As the chemical potential \i increases from 
zero, the critical temperature T c for the chiral symmetry restoration decreases. This critical 
temperature becomes zero at some critical value of the chemical potential fi c = m where the 
chiral symmetry restoration is of first order. One can draw a physical picture of the chiral 
symmetry restoration in terms of overlapping composites. When the temperature is increased, 
chiral condensates begin to overlap as their radius grow up. At some critical point the system is 
mainly composed of overlapping condensates which, as a result, are no longer the good basis for 
describing the thermodynamic properties of the system and the fermionic constituents must be 
taken into account. It is then reasonable to expect that by increasing the baryon density, which 
amounts to increasing the density of chiral condensates, a lower critical temperature would be 
needed for the system to reach the critical point above. In Fig.l we plot the critical lines in the 
(A)-// plane for various values of T. In Fig.2 we plot (A) vs T for various values of ji. In Fig. 3 
we plot the fermion density B/N vs T for various values of \x. For [i > the fermion density is 
discontinuous at a critical temperature T c (/x). 
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4 Imaginary chemical potential 



4.1 The connection with an anyon-like system (anyon confinement /deconfinement) 

Consider the interaction of the Gross-Neveu model (1) with an external U (1) gauge potential 
A^(x). The Euclidean Lagrangian density reads 

£ = -^-ie4)^ i -|(^V < ) 2 . (20) 

where e is the electric charge. Let us consider a constant potential Aq = 9 along the "time" 
direction. We can imagine embedding the model above into a 4-dimensional space. Due to 
the finite length of the xo dimension and the antiperiodic boundary conditions of the fermions 
along it, the system may be viewed as existing in a 3-dimensional hyper-cylinder whose axis 
is the 4 th (unobservable) dimension. The constant potential Aq may now be regarded as the 
"vector" potential generated by a thin solenoid of magnetic flux $ = 9L along the axis of the 
hyper-cylinder. Such a picture corresponds to fermions encircling a thin solenoidal magnetic 
flux and one might expect to encounter Aharonov-Bohm type phenomena [14]. 

The potential in (20) may be gauged away by the transformation 

i) -► e ie9r ip . (21) 

Such a transformation, however, "twists" the antiperiodic boundary conditions for the fermions 
unless 

9 = -Jjfe, k = 0±1,±2,... . (22) 
eL 

The configurations with "twisted" boundary conditions may be viewed as anyon-like excitations 
[14]. For instance, the "quasi-particle" propagator of the above model (20) 

S(x , x; A , (A)) = (p + ie 7o A) + (X))S (x , x; A , (A)) , (23) 

has the following representation as a sum over non-trivial topological paths 

oo 

S (x , x; A ,(X))= ]T e M e^ lL+x °> A °U(x + IL, x; (A)) , (24) 

l=— oo 

where H(xq,x; (A)) are free boson propagators of mass (A). Furthermore, since in our generalized 
representation of the partition function (11) we imposed a constraint on the average fermion 
number, we can have local charge density fluctuations b(r, x) which give rise to the so-called 
statistical gauge field a M (r, x) as 

6(r,x) = ^t( TjX )^( TjX ) _ B = dia 2 (r,x) - d 2 ai(r,x) = e 0ij d iaj (T, x) . (25) 

This (non-dynamical) field would be responsible for the anyon dynamics [22]. Finally, the 
partition function of the free theory i.e (20) at g = 0, has been shown to reproduce the standard 
anyon virial coefficients [23]. 
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The existence of the anyonic excitations above may be tied to a discrete Z-symmetry of 
the 2+1 dimensional Gross-Neveu model interacting with a standard £7(1) gauge field A )l (x) at 
temperature T. The partition function in this case is 

Z = J (PA M )(PV')(^)e~^^ F ^ F ^"^ ( ^" i ^ )V ' l " fW " >l)2] , (26) 

where F^ v = d^A u ] as usual. This theory is invariant under 

A„(t,x) ^ M (r,x)+a M x(r,x), (27) 

^(r,x) e ie ^ T ' x ^(r,x), (28) 

^(r,x) ^(r,x)e- iex(r ' x) , (29) 

which are the usual small gauge transformations provided they are periodic in Euclidean time r 

X (0, x) = X (L, x) and £^(0, x) = d^L, x) . (30) 

In addition to these, the theory is also invariant under large gauge transformations which are 
given by 

X(0,x) = x (L,x) + — k, fc = ±l,±2,.... (31) 

and represent a global Z-symmetry. This symmetry implies the existence of infinitely many 
equivalent Z- vacua in the theory (26). Namely, the gauge field configurations 

AJ(t,x) = (— fc,0,0), fc = 0,±l,±2,..., (32) 

are lowest energy and equivalent since they are connected by gauge transformations of the 
form x(r, x) = (27r/eL)r. Still, such Z-transformations preserve the anti-periodic boundary 
conditions. We therefore conclude that Eq.(32) represents equivalent vacuum configurations of 
(26). 

Going back to (20) provided condition (22) holds, the theory can be interpreted as a Gross- 
Neveu model with a U{\) gauge coupling at finite temperature and lying in one of the Z '-vacua. 
In this case there exist no anyon-like excitations in the spectrum. If condition (22) is relaxed, 
then (20) corresponds to the gauge theory above where the gauge field fluctuates around the 
Z-vacua. In the latter case there exist anyon-like excitations in the spectrum. 

We can now construct an order parameter to distinguish the two cases as follows. Consider 
translational invariant fluctuations around the Z-vacua in (26) such that 

A li = Ajl + 0(T), 9(Q) = e(L). (33) 

Take now the quantity 

(P e -) = (e ie "/o Ld ^W), (34) 

where e/e/ integer and the average is taken with respect to the partition function (26). It 
is not difficult to see that (34) is invariant under the strictly periodic gauge transformations 
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(30). It is not invariant, however, under the Z-transformations (31). To see this consider the 
transformations 

2ir 

6{t) 9{t) + —k , fc = ±l,±2,..., (35) 
eL 

which are symmetries of the theory (26) provided there is Z-invariance. Under (35) 

{Ve)^j 2 * k HVe)^{Ve). (36) 

The physical interpretation of (34) is that Fg = —{1/(5) \a(Vg) gives the free-energy necessary to 
immerse a configuration of charge e in the spectrum of the system [3]. If the theory (26) is in a 
Z-vacuum, then the spectrum can contain only an integer number of fermions or antifermions 
since, as explained earlier, it does not contain anyon-like excitations. In this case one clearly 
expects that Fg = oo or (Vg) = 0. If on the other hand the theory is not in a Z-vacuum, then 
there exist anyon-like configurations in the spectrum and one in general expects that (Vg) ^ 0. 
We conclude that, in this sense, the Z-symmetry of the full gauge theory (26) is connected to 
the "confinement" or "screening" of anyons, where (Ve) is the relevant order parameter. The 
anyonic theory preserves Z-symmetry in the confined phase where (Ve) = and violates the 
Z-symmetry in the deconfined phase where (Vg) ^ 0. 

We are now in a position to give a physical interpretation to the partition function (11) with 
real (9) (imaginary chemical potential), as a result of the following identity 11 

Z{T,V,B) = e-f> F <T> v > B ) = (V S/N ). (37) 

Namely, the free-energy of the Gross-Neveu model with imaginary chemical potential represents 
the free-energy of anyon-like configurations emerging in the spectrum when the model is coupled 
to a £7(1) gauge field which fluctuates around the Z- vacua. 

The above discussion concerning the order parameter (34) parallels the discussion in [13] of 
a similar order parameter for the confinement or screening of incommensurate charges in parity 
invariant QED3 (see also [24]). One of the new observations here is the interpretation of the 
incommensurate charged configurations as being anyon-like. Another important point is that 
the order parameter (Vg) is in fact imaginary for any real value of e. This is also seen from (17) 
or (37). If (34) is to be interpreted as a physical free-energy density, it is necessary that e is 
imaginary. The fact that we obtain an imaginary eigenvalue B = e for the Hermitian operator 
N means essentially that we are not using a positive definite density matrix to describe the 
fluctuations around the Z- vacua of the theory (26). Nevertheless, even in that case we will be 
able to extract useful results for the critical properties of the theory, assuming that the latter 
are universal. 

11 Here we use e = 1. 
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4.2 Chiral Symmetry and Z-vacua 



The chiral symmetry restoration for real (0) can be inferred from the gap equation (15). The 
critical line separating the chirally symmetric from the chirally broken phase in the (0)-T plane 
is obtained by setting (A) = in (15) as 



(0) = — arccos 
Li 



Lm 
1 



(38) 



and is depicted for various values of m in Fig.4 when (9)/T 6 [0, tt]. 



In order to discuss the chiral properties of the theory we shall henceforth consider, for 
concreteness, the massless case where the system is critical already at T = 0. In this case, when 
the chemical potential is real one does not expect any phase transition as the temperature rises 
up. To put it differently and in a more general ground, for real chemical potential the system is 
chirally symmetric at some high temperature and as it cools down chiral symmetry is, in general, 
broken at some lower critical temperature. Having chosen m = this critical temperature is 
T = 0. Consider now the case when (0) is real, where this corresponds to an imaginary chemical 
potential or some fluctuation around the Z-vacua as explained above. Again, at high enough 
temperatures we expect the system to be chirally symmetric. However, chiral symmetry is now 
broken at a non-zero temperature T = 3(0) /2ir. The important point is that as the system cools 
down, chiral symmetry remains broken in a "temperature window" 3(0) /2ir >T> 3(6) /Att and 
then it is restored again. Furthermore, as T — ► 0, it passes through infinitely many "windows" 
of the form ^ 3(0) I 

~2^T l + 3k ~ T - ~2^~ 2 + 3k ' fc -°' 1 ' 2 '-"' ( 39 ) 
in which chiral symmetry is always broken. This essentially means that for any real non-zero 
(0) the T = theory becomes unstable. The relevant physical picture can be read off from Fig. 5 
where we plot (\)/T vs (0) /T in < m/T < 2 In 2. Note that for m/T > 2 In 2 the theory is 
always in a chirally broken phase. 

We can draw a physical picture for the chiral symmetry restoration discussed above. The 
magnetic flux (25) associated with local charge fluctuations catalyzes symmetry breaking as it 
stabilizes the chiral condensates against thermal fluctuations inside the "temperature windows" 
(39) 12 . The fact that in the high-temperature weakly coupled regime even a small imaginary 
chemical potential induces chiral symmetry breaking, shows that the former plays the role of 
a strong catalyst of dynamical symmetry breaking, similar to that of a transverse external 
magnetic field [25] or a constant negative curvature [26]. The chiral symmetry restoration 
transition at T = 3(0) /2ir is of the second-order with mean field critical exponents [14, 20] since 
the order parameter (A) is continuous. The order parameter acquires its maximum value for 

12 This picture may be compared to the one of an external magnetic field in finite temperature 2+1 dimensional 
QED recently discussed in [27]. 
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T = (9} /it and then starts to drop until it reaches zero again at T = 3(9} /Air. The same picture 
holds for all chiral restoration "temperature windows" (39). 

It is also interesting to study the behavior of the free-energy density as a function of the 
temperature for non-zero (9). This is depicted in Fig.6 where we plot the free-energy density vs 
(9)/T. For high enough T the system is in the chirally symmetric phase and to leading- iV the 
free-energy density equals that of N massless free four-component Dirac fermions e.g. Eq.(8). 
At T = 3(9} /2ir the free-energy jumps discontinuously to a local minimum and then raises until 
it reaches a local (positive) maximum at T = 3(9} /An where it drops again discontinuously to 
the value (8). Similar fluctuations occur for all the infinitely many temperature "windows" (39) 



It is rather intriguing that we are able to give analytic expressions for the free-energy density 
of the system at the end-points of the "temperature windows" of chiral symmetry breaking, 
3(9} /2tt, 3(9} /4-7T (in fact this is possible for the end-point of all the "windows" (39)). After 
some algebra we obtain respectively 



where, Cl2(<fi) = ImfL^e 1 ^)] is the Clausen function [21]. Note that Cl2(it/3) is the absolute 
maximum of this function which is a well-documented irrational number [28]. 

The results for the free-energy density have been on purpose written as above, in order to be 
compared with the expected scaling form of the free-energy density of a conformal field theory 
(CFT). The point is that the chiral phase transition in the 2+1 dimensional model above is of 
second-order and one would expect it to correspond to the universality class of a 3-dimensional 
CFT. When a CFT is put in a finite-size geometry (i.e. in a slab with one finite dimension of 
length L), its free-energy density scales as [18] 



The quantity c coincides with the central charge in d = 2 and has been recently proposed [29] to 
be a possible generalization of a C7-function in d > 2. Clearly, from (40), (41) and (42) we see 
that the chiral transition above appears to be connected to new 3-dimensional CFTs. The fact 
that c in (40) is less than the corresponding free-field theory value (8) and it is negative in (41), 
seems to imply that the above critical theories may not be unitary. Nevertheless, such theories 
may conceivably correspond to three-dimensional versions of the non-unitary two-dimensional 
Lee- Yang model [32]. 

The middle point of the chiral transition "temperature window" (39) is also interesting. At 



as T -► 0. 




(40) 



(41) 
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this point (A) takes its maximal value. From (15) we obtain for (9) = ir/L, 

This value equals the value one obtains for the mass of the elementary bosonic modes in the 
2+1 dimensional O(N) vector model at finite temperature, when the theory is critical at T = 
[30]. Plugging this into (17) we can calculate the free-energy density which, by virtue of some 
non-trivial polylogarithmic identities [30, 31] is found to be 

h = 16 CR 

N 5 2ttL 3 y ' 

This is exactly minus the free-energy density of the 2+1 dimensional 0(4N) vector model at its 
non-trivial critical point. In principle, a positive free-energy density which would correspond to 
negative pressure and negative entropy seems to be a rather unphysical result [33]. However, 
one might try to construct a physical system where the result (44) could make sense. This is the 
0(AN) M = 1 supersymmetric sigma model [34, 35] in the presence of an external Aq potential 
at finite temperature. To leading- N, the free-energy density of this theory is simply given by the 
sum of the free-energies of fermions and bosons. In general, supersymmetry is expected to be 
broken at any finite temperature [36]. Nevertheless, it may happen that for some temperature 
the fermion contribution, given by (44), and the boson contribution cancel each other and the 
system becomes supersymmetric again. Note that the matching of the bosonic and fermionic 
degrees of freedom is correct - the N = 1 supermultiplet in three-dimensions requires 2N two- 
component Majorana fermions. A similar bosonization of fermions has been recently discussed 
in [37]. 



5 Summary and Discussion 



In this work we studied the 2+1 dimensional Gross-Neveu model in the presence of an 
imaginary chemical potential and argued that it provides an interesting ground for probing the 
properties of chiral symmetry at non-zero temperature. In Sec. 3 we proposed a generalization of 
the well-known formula for the canonical partition function and presented analytic expressions 
for the free-energy density, the fermion number density and the chiral order parameter (A). We 
demonstrated that this general formalism includes the standard grand-canonical formalism when 
the chemical potential is real. In Sec. 4 we focused on the case of imaginary chemical potential. 
We considered a U{\) gauge field coupled to the Gross-Neveu model at finite temperature. Such 
a theory possesses infinitely many equivalent Z-vacua. We showed that when the gauge field 
fluctuates around these Z-vacua, as given by Eq.(33), anyon-like excitations manifest themselves 
in the spectrum of the 2+1 dimensional Gross-Neveu model. This is established by showing that 
the free-energy necessary to immerse such an excitation is non-zero. Furthermore, we provided 
evidence that the expectation value of the abelian Polyakov loop, which is the order parameter 
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for establishing whether the above U(l) gauge theory resides in one of the Z-vacua, coincides 
with the canonical partition function for imaginary values of the fermion number density. This 
way we gave a physical interpretation to the canonical partition function with imaginary chemical 
potential as being the free-energy needed to immerse an excitation of imaginary charge in the 
spectrum of the 2+1 dimensional Gross-Neveu model. The latter excitation is related to the 
anyon-like excitations discussed above. Finally, we gave some results connected to the properties 
of chiral symmetry in the presence of an imaginary chemical potential at finite temperature. 

Studies of the critical thermodynamic properties of a system undergoing a phase transition 
are intimately connected to the Lee- Yang zeroes [38, 39]. These are the zeroes of the partition 
function for imaginary values of the external magnetic field, the latter being the "conjugate" 
variable of the relevant order parameter. It is then conceivable that one could try to investigate 
the critical properties of a system at finite density by studying its partition function (or its 
free-energy density), at complex values of the number density, the latter being the "conjugate" 
variable of the chemical potential. In this sense, the results presented in this work are closely 
related to a Lee- Yang zeroes analysis. 

These results can be extended in several directions. For example, it is possible in principle 
to study numerically the partition function (11) for complex chemical potential and establish 
the chiral symmetry restoration "windows" advocated above. It is also possible to numerically 
compute the free-energy densities at the critical points and compare them with our analytic 
results (40) and (41). Furthermore, it would be interesting to compare the free energy density 
(17) to the one recently proposed by Laughlin [40] in order to reproduce the anomalous behaviour 
of the thermal conductivity in some high-temperature superconductors of the BSCCO family [41] 
13 . Since these compounds exhibit a new kind of phase transition induced by strong catalysts of 
dynamical symmetry breaking, it is plausible that our generalized canonical partition function 
formalism can be used to explain such unusual behavior. Finally, based on the results of the 
present work, it may be interesting to explore the possibility of finite temperature supersymmetry 
restoration in = 1 supersymmetric theories [43]. 
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Figure 1: Critical lines in the (A) 
such that T c « 0.72. 



— ijl plane for various T < T c . 
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For simplicity we chose 



OC 




T 




Figure 4: Critical lines in the (6) — T plane for various m when (9)/T € [0,7r]. The line 
(6) = (2tt/3)T is an asymptote. 
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Figure 5: (A)/T vs (9)/T for < m/T < 2 In 2. For fixed, non-zero (0), T decreases along the 
horizontal axis as one moves to the right. 
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Figure 6: The free-energy density vs (9)/T. As in Fig.5, T decreases along the horizontal axis 
as one moves to the right. 
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